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FORMAL  I SM/THEORY 


BACKGROUND 

Following  Tatarskii1,  we  consider  the  angular  jitter  measured  by 
a  lens  type  antenna  where  the  scattered  fields  are  given  by  the  Rytov 
approximation.  The  use  here  of  the  Rytov  approximation,  which  is  a  form 
of  weak  scatter  theory,  is  justified  in  Appendix  A.  There  it  is  shown 
that  for  the  propagation  paths  considered,  the  estimated  amplitude  fluctu¬ 
ations  are  small  compared  to  the  level  at  which  weak  fluctuation  theory 
breaks  down.  The  development  of  a  general  equation  capable  of  predicting 
results  for  the  variation  of  all  parameters  of  interest  will  be  described. 
A  Fortran  code  was  written  and  validated  to  obtain  results  from  this 
general  formula. 

The  modeling  geometry  is  shown  in  Figure  1.  The  radar  antenna 
is  immersed  in  a  turbulent  layer  of  dust  and  is  tracking  a  target  in 
angle.  The  angle  estimator  is  assumed  to  correspond  to  the  angular  posi¬ 
tion  of  the  target  image  intensity  pattern  in  the  image  plane  of  the 
antenna.  Due  to  the  fluctuations  induced  by  the  turbulent  propagation 
medium,  the  position  of  the  target  image  will  fluctuate.  Tatarskii1  has 
considered  this  problem  for  the  case  of  an  aperture  with  uniform  illumina¬ 
tion.  We  have  modified  his  derivation  to  include  the  effect  of  an  aperture 
amplitude  taper. 

From  Appendix  B,  the  mean  squared  deviation  in  the  angle-of- 
arrival  (in  one  plane  only)  is 
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Propagation 


R  =  Radius  of  antenna  aperture 
6  =  Target  elevation 

P  =  Distance  measured  transverse  to  the  propagation  path 
Z  =  Altitude 

L  =  Path  length  through  the  dust  layer 


Figure  1.  Modeling  geometry. 
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(1) 


2R 

J  dp  p  I(p) 
o 


d2P(p)  +  1_  dD(p)~ 

dp2  p  dp 


where  k  =  free  space  RF  wave  number. 

Ae  and  I ( p )  depend  on  the  aperture  design  taper  or  weighting 
and  are  calculated  from  equations  (B4),  and  (BIO)  respectively.  P(p)  is 
the  phase  structure  function  evaluated  at  the  antenna  aperture  as  a  func¬ 
tion  of  distance  (p)  transverse  to  the  propagation  path. 

POWER  SPECTRAL  DENSITY  FUNCTION 


We  assume  a  power  spectral  density  (PSD)  which  is  locally  homo¬ 
geneous  and  isotropic  and  an  index  of  refraction  variance  which  is  a 
function  of  altitude. 


*n(K,Z)  = 


Lo3  r<N/2>  An2(Z) 

TT3/2  r( N/2-3/2 )  (1+L  2K2)N/2 


=  *J(K)  An2(Z)  (3) 

where 

4>n(K,Z)  =  PSD  of  index  of  refraction  fluctuations 

An2(Z)  =  index  of  refraction  variance  of  the  air-dust  medium  as 
a  function  of  altitude 

L0  =  outer  scale  of  turbulence 

K  =  spatial  wave  number 

N  =  3-D  power  law  exponent. 
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Using  equation  (3)  and  equations  (1.50)  and  (8.11)  from  Refer¬ 
ence  2,  D(p)  may  be  expressed  as 


D(p)  =  4tt2  k2  J  dK  K  [l-J0(Kp)]  (4) 

0 

•  /  dX'  An2(X' )  il+cos[-2(L.lXH]j 
o  (  k  | 

where  J0(Kp)  =  Bessel  function  of  the  first  kind  of  zero  order,  and  the 
integration  in  X'  is  over  the  propagation  path  from  the  layer  boundary  to 
the  antenna  at  slant  range  L  away.  (See  Figure  1) 

Making  a  change  of  variables  from  X'  to  X  such  that  X  =  L-X'  and 

using 


d2J  (Kp) 

0 _  +  _ 

dp2  p 


dJQ(Kp) 

dp 


-  -K2Jq(Kp), 


we  may  differentiate  D(p)  to  obtain 


d.2£ipJ.  +  1  =  4-rr2k 2  /  dK  K3  JQ(Kp)  *°(K)  (5) 

dp2  p  dp  o  1 

.  JL  dX  An2(x)  [1  +  cos  (*!*)]  . 

o  k 

Here  the  integration  in  X  is  from  the  antenna  to  the  layer  boundary. 
Substituting  (5)  into  (1)  we  obtain 

.  3  2R  oo 

aa2  =  —  f  dp  p  1(e)  /  dK  K3  4>°(K)  Jq(Kp)  H(K)  (6) 

e  0  0 
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where 


H(K)  =  J  dX  An2(X)  [1  +  cos(i^)] 
o  k 


LOADING  PROFILES 


(7) 


Let  the  target  elevation  angle  be  e  with  z,  =  csc(e),  and  the 
antenna  height  Zr.  For  uniform  loading  over  the  layer 

H(K)  »  An2  J  dX  [1  +  cos  {—) )  =  An2  L[1  +  _L  sin(i!il)]  (8) 
ou  o  k  ou  K2L  k 

where  L  =  c(Z^-Z^),  and  Z^  is  the  height  of  the  dust  layer. 


For  exponential  loading 
“  -2Z/hs 


H(K)  =  An2e  z  j  dZ  e 
ZR 


Anz 


ch$  e 


-2ZR/hs 


[1  *  Cost'Ll)] 


cov-I:"cZr  )  -  (K-lch.s.)  s  in  ] 

1  +  k  J  1  2k  J  v  k  _ 

'  2k  J 


(9) 


where  the  mean  and  also  the  root-mean-square  deviation  in  the  index  of 
refraction  are  assumed  to  vary  as  exp(-Z/hs)  so  that  Ane2  varies  as 
exp(-2Z/hs).  Anou2  and  A noe2  in  equations  (8)  and  (9)  are,  respec¬ 
tively,  the  index  of  refraction  variances  in  the  dust  laden  air  medium  for 
the  uniform  and  exponential  layers  at  zero  altitude. 


VARIANCE  OF  THE  INDEX  OF  REFRACTION  OF  THE  MEDIUM 


Following  Reference  3,  we  assume  that  the  mean  value  for  the  in¬ 
dex  of  refraction  of  the  medium  is  given  by  the  Clausius-Mossotti  formula 
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(10) 


n  =  l+ni  .  1*1 

2  \nb2+2/  pb 

and  that  the  root-mean-square  deviation  in  n  is  equal  to  1/2  of  n1 

An  =  1/2  nx  .  (11) 

Here  nb  is  the  index  of  refraction  of  the  bulk  material  which  can  be 

closely  approximated  by  using  its  real  part  only,  p  ,  p.  are  the  densi- 

m  D 

ties  of  the  medium  and  bulk  material,  respectively.  The  assumption  An  = 
1/2  nt  corresponds  to  the  worst  case  situation,  which  is  expected  at  early 
times  after  the  pedestal  cloud  has  formed.  As  the  turbulent  energy  in  the 
cloud  is  transferred  into  heat,  the  factor  of  1/2  will  eventually  become 
much  smaller. 


Substituting  (10)  into  (11)  we  obtain 


An  =  3/4 


(12) 


This  may  be  restated  in  a  convenient  form  by  making  use  of  the  fact  that 
pm  is  generated  by  a  scoured  layer  of  thickness  aZ  of  the  bulk  material 
whose  density  is  p^.  Then  for  uniform  loading 


An 


ou 


3 

4 


and  for  exponential  loading 


An 


oe 


(13) 


(14) 


The  variance  of  the  index  of  refraction  of  the  medium  is  then  obtained  by 
taking  the  square  of  (13)  or  (14)  as  desired. 

WORKING  EQUATION 

The  parametric  investigation  was  performed  through  the  evalua¬ 
tion  of  equation  (6)  which  we  restate  by  substituting  for  <&°  from  equation 
(2) 

2*3/2  L  3  r(N/2)  2R  - 

a  2  =  - - -  x  /  dp  p  I(p)  J  dK  K3  J  (Kp)  H(K)  (1+L  2K2)‘N/2 

“  A  2  r(N/2  -  3/2)  o  o  °  0 

(15) 

where  Ae  and  I(p)  are  given  by  equations  (B12)  through  (B17)  and  H(K)  is 
taken  from  (8)  or  (9)  as  desired.  The  two  integrations  in  equation  (15) 
were  computed  numerically. 
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RESULTS 


A  list  of  parameters  affecting  the  angle-of-arrival  problem 

which  have  been  studied  is  given  below.  Note  there  are  two  parameters 
(antenna  aperture  illumination,  radar  height)  which  the  designer  could 

possibly  control.  (Changing  the  environment  by  relocating  the  radar  site 

or  by  modifying  the  terrain  from  which  the  dust  is  lofted  are  other 
obvious  considerations.) 

.  CHARACTERIZATION  OF  THE  ENVIRONMENT 

-  PSD 

Turbulence  Outer  Scale  Size  (L0) 

3-D  Power  Law  Exponent  (N) 

Bulk  Index  of  Refraction  (n^) 

-  Dust  Loading  Profile 

Uniform 

Exponential 

Scoured  Layer  Removed  (aZ) 


.  RADAR  OPERATION 

-  Antenna  Aperture  Effects  (Illumination  Functions) 

-  Target  Elevation  (e).  Radar  Height  (Zr) 
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NOMINAL  CASE 


Unless  otherwise  specified,  we  assume  the  following  nominal 
conditions: 


Exponential  loading,  hs  =  6  m 
2  cm  of  caliche  scoured  (n^  =  2.2) 
Kolmogorov  spectrum,  N  =  11/3 
L0  =  10  m,  e  =  90* 

Uniform  aperture  illumination 
Radar  height,  Zr  =  0  m 

For  these  conditions  we  find  oa  «  1.7  mrad. 

VARIATIONS  ABOUT  THE  NOMINAL 


Figure  2  demonstrates  the  variation  in  aa  with  outer  scale  and 
the  3-D  power  law  exponent  of  the  PSD  for  an  exponentially  loaded  layer. 
Variations  in  attributable  to  uncertainty  in  Lq  can  approach  factors  of 
±  2  for  l<Lo<100  m  for  a  given  N.  For  a  fixed  choice  of  outer  scale,  a 
variation  in  N  from  3  to  5  can  produce  variations  of  ±  /3  in  aa. 


Figure  3  is  a  repeat  of  Figure  2  except  that  the  2  cm  of  scoured 

caliche  has  been  distributed  uniformly  over  a  100  m  layer.  It  is  observed 

that  the  values  of  oa  predicted  for  the  uniform  layer  are  about  a  factor 

of  3  smaller  than  for  the  exponential  layer  with  a  6  m  scale  height.  This 

can  be  explained  as  follows.  Let  the  mean  turbule  diameter  for  a  qiven  L 

3  o 

and  N  be  equal  to  d.  Then  in  traversing  a  layer  thickness  Z,  aa2  will  be 
proportional  to  the  variance  in  the  index  of  refraction  for  an  average 
turbule  (An2)  times  the  number  of  turbules  passed  (m).  Thus 


a,u 


a,e 


An 


An. 


(16) 
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a.(MRAD) 


a.(MRAD) 


where  the  subscripts  u  and  e  refer  to  uniform  and  exponential  respec¬ 
tively,  For  a  constant  layer  «  Z^/d  and  for  an  exponential  layer  * 
hs/d.  Also*  referring  to  equations  (13)  and  (14) 


An 


2 

u 


An 


e 


2 


Then 


a, U 


a,e 


VZh 

zl\ 


1 

e  _ 


4 


(17) 


(18) 


approximating  the  numerical  result  of  1/3  quoted  earlier.  The  results  for 
the  uniform  layer  have  been  compared  with  a  previous  result  published  by 
Thompson  in  Reference  4.  They  are  found  to  be  in  excellent  agreement  when 
the  differences  in  the  scoured  layers  assumed  are  taken  into  account. 


In  Figure  4  we  have  investigated  variations  of  oa  with  scale 
height  for  exponential  loading.  The  increase  in  aa  with  decreasing 
scale  height  may  be  explained  by  the  same  reasoning  as  was  used  above  to 
compare  the  results  for  the  uniform  and  exponential  layers.  The  variation 
of  hs  about  the  nominal  6  m  value  results  in  fluctuations  of  a  factor  of 
*  0.5  to  1.7  in  oa. 

Figure  5  shows  the  effects  of  target  elevation  and  also  radar 
height  on  oa.  For  a  decrease  in  target  elevation  from  the  nominal  case 
of  90°  down  to  15"  we  see  a  factor  of  about  2  increase  in  oa.  This  rise 
results  from  an  increase  in  path  length  through  the  dust  with  decreasing 
look  angle.  For  a  fixed  target  elevation  and  our  nominal  6  m  dust  scale 
height,  we  note  a  decrease  in  aa  by  a  factor  of  «  .7  when  the  antenna 
face  is  lifted  2  m  above  the  ground.  This  decrease  in  o  is  caused  simply 
by  the  fact  that  a  large  quantity  of  the  former  propagation  environment  is 
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a,(URAD) 


Figure  4.  oa  vs  dust  density  scale  height. 
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o.(MRAD) 


TARGET  ELEVA TION(DEG) 


Figure  5.  aa  vs  target  elevation  and  radar  height. 
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now  below  the  antenna  and  has  no  effect.  It  should  be  realized  that  for 
smaller  scale  heights  this  factor  would  become  more  significant,  and  for 
large  scale  heights  the  effect  of  a  2  m  radar  height  would  be  negligible. 


The  three  separate  aperture  amplitude  design  tapers  are  plotted 
in  Figure  6,  with  their  corresponding  peak  sidelobe  levels  also  listed  for 
reference.  The  resultant  variations  in  aa  as  a  function  of  design  taper 
and  outer  scale  are  shown  in  Figure  7.  It  is  seen  that  with  increasing 
amplitude  taper,  corresponding  to  greater  sidelobe  suppression,  there 
results  an  increase  in  aa.  Thus  there  exists  a  tradeoff  between 
sidelobe  suppression  and  angle  of  arrival  jitter.  The  peak  increase  in 
aa  is  a  factor  of  «  1.6.  This  increase  occurs  because  a  tapered 
aperture  gives  lower  sidelobes  at  the  expense  of  a  somewhat  broader  beam 
and  a  broader  beam  in  space  lets  in  rays  from  a  wider  range  of  angles. 
The  possible  significance  of  this  effect  remains  to  be  investigated. 


The  effect  of  varying  material  properties  when  the  scoured  layer 
thickness  is  held  constant  is  obtained  directly  from  equation  (13)  or  (14) 


oa  *»  An  » 


%2-i 

n"b2+2 


(19) 


Table  1  lists  the  bulk  indices  of  refraction  and  resulting  variations  in 
(nj^2 -1 ) / ( n|32+2 )  for  representative  materials  of  interest.  We  note  a 
resultant  factor  of  «  .6  to  1.4  in  aa  about  the  nominal  case. 


Table  1.  Material  Variations.  (From  Gutsche's  data.  Reference  5.) 


Material 

nb 

nb2-l 

nbz+2 

Dry,  sandy  soil 

1.58 

.33 

*  Caliche 

2.20 

.56 

Wet  Clay 

3.34 

.77 

*  Nominal  Case 
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B  T" - —  . — - 


p/R 


Sidelobe 
Level  (dB) 

(T)  F(p)  =  Const  18 

©  F(p)  =  (1-P2/R2)  25 

©  F(p)  =  (l-p2/R2)2  31 


Figure  6.  Aperture  amplitude  design  tapers. 
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OUTER  SCALE(U) 


(T)  F(p)  =  Const 
(D  F(p)  =  (l-p2/R2) 
(D  F(p)  =  (1-p2/R2) 2 


Figure  7.  oa  vs  aperture  taper  design  and  outer  scale. 


CONCLUSIONS 


Results  are  summarized  in  Table  2.  It  is  seen  that  for  the 
nominal  conditions  of  2  cm  of  scoured  surface  we  predict  about  2  mrad  of 
jitter.  However,  notice  that  aa  is  directly  proportional  to  the  amount 
of  surface  removed  and  this  quantity  is  highly  uncertain.  Variations  in 
other  quantities  lead  to  factors  of  -  ±2  about  the  nominal  case.  In  addi¬ 
tion,  with  specially  selected  (worst  case)  parameter  variations,  as  much 
as  10  mrad  of  jitter  is  possible  with  the  2  cm  of  scoured  surface. 

It  is  appropriate  to  consider  the  amount  of  signal  attenuation 
predicted  for  this  dust  loading.  From  Reference  5,  Figure  1  we  note  that 
for  2  cm  of  surface  removed,  there  is  less  than  15  dB  of  2-way  signal 
attenuation  for  the  worst  case  combinations  of  material  composition  and 
its  size  distribution  in  the  medium.  Thus,  a  couple  of  milliradians  of 
jitter  may  be  possible  under  conditions  for  which  attenuation  is  not 


severe. 


Table  2.  Results. 


VARIATIONS  RELATIVE  TO  NOMINAL  CASE 

o  (nom)  «  1.7 

a 

mrad 

aa 

(min) 

aa  (max) 

QUANTITY  VARIED  aa 

(nom) 

aa  (nom) 

PSD 

.3 

1.4 

Loading  (scoured  layer  removed) 

Directly  Proportional 

Material  properties 

.6 

1.4 

Scale  Height 

.5 

1.7 

Radar  Height 

.7 

1.0 

Target  Elevation  1 

.0 

2.0 

TRADEOFF  BETWEEN  SIDELOBE  SUPPRESSION  AND  JITTER 
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APPENDIX  A 

JUSTIFICATION  FOR  THE  USE  OF  WEAK  SCATTER  THEORY 


The  formalism  which  has  been  developed  utilizes  results  from 
weak  scatter  theory.  According  to  Reference  6  weak  fluctuation  theory  is 
valid  when  the  log-amplitude  variance  is  less  than  about  .2  to  .5.  For 
the  case  of  a  Kolmogorov  spectrum,  the  log-amplitude  variance  is  given  by 


2  _ 


=  .586 


onWl/6 

_  i  2/3 


(Al) 


for  uniform  loading  and  a  path  length  L.  Vie  approximate  the  "ominal 
exponential  layer  (hs  =  6  m)  by  a  uniform  layer  of  height  m.  t, **■'  is 
taken  to  correspond  to  that  at  ground  level  for  the  nominal  conditions. 
This  should  yield  a  conservative  estimate.  Then 


L  =  6  m 

L0  =  10  m 

k  =  (2ir/.03)m“ 1 

An  -  1  lHl±) 

4  \  2.2z+2  /  6 

and 

ox2  -  3.4  10-3 

This  is  two  orders  of  magnitude  smaller  than  the  value  for  which  weak 
fluctuation  theory  breaks  down.  It  is  safe  to  assume  that  ox2  will  be 
smaller  than  .2  for  all  cases  we  have  considered,  so  that  the  use  of  weak 
scatter  theory  is  justified. 


BLANK  PAGE 
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APPENDIX  B 


ANGLE  OF  ARRIVAL  FLUCTUATIONS  FOR  AN  ANTENNA 
WITH  A  TAPERED  APERTURE  AMPLITUDE  DESIGN 


As  stated  earlier,  we  have  assumed  that  the  relevant  angle 
estimator  is  that  of  the  angular  coordinates  of  the  "center  of  gravity"  of 
the  intensity  pattern  of  the  target  image  in  the  focal  plane  of  the 
antenna.  Following  Tatarskii1,  the  angular  deflection  of  the  center  of 
gravity  of  the  image  (in  one  plane  only)  is  given  by 

o0  =  I  //  Im  [*(n,,C)  !iLals£l]  F2(n' »£' )d€'dn' 

_ _ I - in! _  (bi) 

//  )|  2  F2(n',c)  dn'dC 


where 

k  =  free  space  RF  wave  number 

i|i(n'»£')  =  scattered  field  at  the  antenna  aperture  in  the  Rytov 
approximation 

F(n',c')  =  antenna  aperture  amplitude  design  taper 
Z  =  aperture  area 

(This  may  be  obtained  from  Reference  1  page  287,  equation  (13)  by  letting 
ip-*-F  ip . ) 


Letting  =  Aq  exp(x+iS),  where 

x  =  log-amplitude  of  the  scattered  field 
S  =  phase  of  the  scattered  field 

we  may  write  equation  (Bl)  as 
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//  F2(n\C')exp(2X)  ?S(  n.'.» il^n'dg 
-1  Z  3n' 

k  //  F2( n' , 5‘  )exp(2X)  dn'dc'  ~ 

z 


( B2) 


It  is  seen  in  equation  (B2)  that  the  main  contribution  to  oq 
is  from  phase  fluctuations  since  a0  =  0  for  S  =  constant.  Then  to  first 
order  we  may  neglect  the  effect  of  amplitude  fluctuations  and 

ct .=  —  !/  F2(n\5’)  dn'dC  (B3) 

0  kAg  z  3n' 


where 


Ae  =  "effective  area" 

=  ff  F^n'.C’Jdn'dc' 

z 


(B4) 


The  mean  squared  deviation  in  a  is  then 


Oa2  =  <«02> 


-L-  ff  jj  F2(n',C'  )F2(  n",  £")<  —  n-?-C< )  1S.— \  > 

k2A  2  z  z  3n'  3n" 

e 


dn'd^'dn"dc" 


( B5) 


The  expectation  of  the  phase  derivatives  may  be  written  in  terms 
of  the  phase  structure  function  as 


,3S(n\C‘)  8S(n'\C") 


>  = 


3n‘ 


3n" 


<S( n' * C1 )  S(n",£")> 


3n'  3n" 


1 


2  2 

— - - B(n'-n\  C-V')  = 

8n'9n"  3n' 2 

1  _if_  D(n'-n,,,C,-C") 

2  an'2 


(B6) 


where  it  is  assumed  that  the  phase  correlation  function  B  and  the  phase 
structure  function  D  depend  only  on  the  distance  between  the  points 
(n',£‘)  and  (n",S").  Substituting  (B6)  into  (B5),  we  may  write 


2  =  - —  If  ff  F2(n',C')  F2(  n",  £" ) 


a  2k2A  2  Z  Z 
e 


an' 


dn,dc'dn"dc" 


(B7) 


We  now  make  a  change  of  variable  from  (ri"»C")  to  (n,£)  such  that  n  = 
n'-n"»  S  =  C'-£".  Keeping  in  mind  that  the  area  z  is  a  circular  aperture 


of  radius  R, 

1___  ff 

dndf;  a2P(n»c) 

a  2k2Ag2  /  n2+C2  <2R 

an2 

.  ff  dn'dc1  F2(tV,£')  F2(V-n,C'-0  •  (B8) 

z 

Changing  to  polar  coordinates  (p,4>)  and  assuming  that  the  structure  func¬ 
tion  is  locally  isotropic  so  that 

D(n,0  *  P  (  /n2+S2)  and 

a2P(p)  _  d2P( p)  +  ji2  dP( p) 
an2  p2  dp2  p^  dp 

-  c«S**2f£t£l  .  1  «<£» 

dp2  p  dp 
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we  obtain 


on  p 

o  2  = _ !: _ /  dp  p  /  d<p  [cos2<p  £  +  2.  sin2$  2Pi_p-2]  I ( p)  (B9) 

°  2k2Ag2  o  o  dp2  p  dp 

where 

I ( p )  =  "illumination  function" 

cos_1(p/2R)  R  , - 

=  4  /  d *'  /  dp'  p1  F2(  p' )  F2[  /(p*- p)2  ]  -  (BIO) 

o  p/2cose 

Performing  the  integration  over  $  in  equation  (B9)  we  have 

2R 

0_2  =  -  JL ,  /  dp  P  I(p)  [d  D(p)  +  I  dP(p)-  ]  .  (811) 

p  dp 

of  equation  (BIO)  is  shown  in 
(BIO)  for  3  separate  amplitude 

tapers* 


2k‘A  * 
e 


The  geometry  for  the  integration 
Figure  Bl.  We  evaluate  equations  (B4)  and 


Case  1: 

F(p)  =  1 

Ag  =  nR2  (B12) 

I(p)  =  2R2  [cos-Mp/ZR)  -  ( p/2R)  /  1-p2/4R2]  (B13) 


*  I(p)  was  evaluated  analytically  for  Case  1,  and  numerically  for  Cases 
2  and  3.  For  Cases  2  and  3  I(p)  was  then  expressed  in  terms  of  poly¬ 
nomial  expansions  through  use  of  a  curve  fitting  routine. 
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Case  2: 


R(p)  =  l-p2/R2 

Ag  =  uR2/3  ( B14) 

I ( p )  =  R2  [A2+B2(p/R)  +  C2(p2/R2)  +  D2( p3/R3) 

+  E2(pVR4)  +  F2(p5/R5)  ]  ( B15) 

Case  3: 

F(p)  =  (l-p2/R2) 2 

A  =  ttR2/5 

e 

I(P)  =  R2  [A3+B3(p/R  +  C3( p2/R2)  +  D3(p3/R3) 

+  E3( p4/R4)  +  F 3( p5/R5)  ] 

where 

A2  =  .6220536 

B2  =  .1518489 

C2  =  -1.859473 
D2  =  1.699667 

E2  =  -.5666203 

F2  =  .06177243 


(B16) 

(B17) 
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